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The upper critical field problem in MgB2. 
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The upper critical field Hc2{T) for MgB2 is analyzed in terms of single and multi-band Eliashberg 
models. The relatively high value of Hc2{0) ~ 14 to 18 T can be understood, if a strongly coupled 
subgroup of heavy quasiparticles is involved in the superconductivity. The sizable coupling to both 
low- and high-frequency bosons is essential. This picture is supported by LDA calculations of Fermi 
velocity distribution over different sheets of the Fermi surface with special emphasis on the hole 
tubes. The possible origin of anharmonic soft modes is briefly discussed. Similarities and differences 
with transition metal borocarbides are figured out. 



The recent discovery of superconductivity in MgB2 ||] 
has initiated an immediate broad research activity due to 
the high transition temperature Tc ~ 40 K in a seemingly 
ordinary s — p metal. The present discussions mainly fo- 
cus on measured quantities Tc and the gap 2Ao with re- 
spect to the topics: (1) weak or strong coupling in terms 
of standard BCS theory, (2) its possible relationship to 
the superconductivity in transition metal borocarbides, 
(3) symmetry of the order parameter, and (4) electron- 
phonon (el-ph) interaction or Coulomb repulsion based 
mechanisms of pairing In context of the stan- 

dard phonon mechanism: (i) a dominant role of inter- 
mediately coupled high-frequency boron phonons is as- 
sumed (the simple "metallic hydrogen scenario"); or (ii) 
a strong coupling scenario is discussed. In all those dis- 
cussions, essentially the validity of the standard isotropic 
single band (ISB) picture is implicitly assumed. How- 
ever, the low temperature value of the upper critical field 
Hc2(0), the fundamental quantity of a type-II supercon- 
ductor, has not analyzed so far. In this context we would 
like to point out an interesting puzzle: with the com- 
puted and commonly accepted average band Fermi veloc- 
ity i;i?=:8.9xl0^ cm/s and the plasma frequency ujpi—7 
eV, the resulting clean-limit London penetration depth 
Al = c/LLipi—29 nm nearly equals the BCS coherence 
length ^0 = 0-l8fivF/kBTc=29 nm while experiments on 
He and Hc2 yield a Ginzburg - Landau parameter k=26 

The ISB model M is the most developed part of the 
theory of superconductivity. It describes quantitatively 
the renormalization of the physical properties of metals 
due to the electron-phonon interaction. The input mate- 
rial parameters are the density of states at Ep, N{0), the 
Fermi velocity vp, the impurity scattering rate ^imp, the 
paramagnetic impurity scattering rate 7™, the Coulomb 
pseudopotential /x*(« 0.1), and the electron-phonon 



spectral function a'^F{fl). The temperature dependen- 
cies of physical properties are determined mainly by the 
first inverse moments of the spectral function - the cou- 
pling constant X—2 f dfla^ F / ^l and the average boson 
energy Qiog = exp[(2/A) dfl log{n)a^F{n)/n]. 

Our analysis is based on numerical calculations of 
i7c2 (r) solving equations first presented by Prohammer 
et al. 1^. In the intermediate and weak coupling cases, 
the clean limit value of the calculated orbital -ffc2(0) prac- 
tically does not depend on the shape of the spectral func- 
tion 1^,0 . In the spirit of Langmann ||^ the s-wave clean 
limit ISB-iJc2(0) can be approximated by the following 
factorized formula 

-ffc2(o)[T] « o.o2r2[K](i + A)2-V4[ioWs]. (1) 

It is convenient to express v^p in Eq. |^ using the re- 
lated experimentally available quantities: the cell vol- 
ume V, the Sommerfeld constant (per mole) 75 = 
7r^fc^A^(0)(l + X)NAV/'i and the bare plasma frequency 
^ WiV(0)4/3 



nNAVkl{l + \)^l 



(2) 



Substituting Eq. | into Eq. |l| one arrives at the relevant 
ratio 



7s[mJ/mol x K2]T2[K](1 + A)i-' 



(3) 



the value of which must be Q «1 in a ISB-like case. A 
strong deviation Q from unity indicates the inapplicabil- 
ity of the s-wave ISB model. 

To proceed within the Q-check, let us estimate the cou- 
pling constant A and the corresponding averaged boson 
frequency ($7). Recent neutron data H show phonon 



density of states (PDOS) extending to above 100 meV 
with main peaks at about 90, 55, 37, and 32 meV, as 
well as an unexpected anomalous soft mode near 17.5 
meV below the maxima of acoustic modes. The measured 
PDOS is in agreement with Q]j=65 meV, estimated from 
specific heat data [0. The values for zone-center opti- 
cal mode obtained from density-functional theory based 
calculations 58 meV Q, 62 meV ||l|, 73 meV ||l| are 
rather similar and for our goal we adopt an Einstein spec- 
trum with ri=60 meV. With the standard value /i*=0.1 
at the cutoff frequency of 6(^2) one needs A=0.86 to ob- 
tain Tc=40 K. If the dominating coupling were that by 
the 17.5 meV feature A=2.5 would be needed. 

Now we are in a position to apply our Q-test to selected 
superconductors. The results are presented in Table 1 |. 
For the weakly anisotropic transition metal Nb it is ap- 
proximately fulfilled, whereas strong deviations do occur 
both for the transition metal borocarbide YNi2B2C and 
for MgB2 under consideration. The latter one is checked 
in the weak and as well as in the strong coupling cases. 

These difficulties can be resolved in the framework 
of a more complex effective iV-band model. For an 
anisotropic system it is obtained by dividing the Fermi 
surface into N parts and approximating the k-dependent 
quantities in each part by their mean values p^ . The 
input material parameters of such a model are the den- 
sities of states at Ep Ni{0), the Fermi velocities vp^i, 
the impurity scattering rates ^imp,ij , the Coulomb pseu- 
dopotentials /^*(« 0.1), and the electron-phonon spectral 
functions a^^Fiuj). For the sake of simplicity, as a first 
step, we consider here such effects in the frame of an 
effective two-band model. Recently it was successfully 
applied to transition metal borocarbides in a combined 
study of Eliashberg theory and de Haas van Alphcn effect 
0. Unfortunately, due to the lack of single crystals, at 
present there is no available experimental information on 
the electronic structure near the Fermi surface of MgB2. 
For this reason, we performed electronic structure calcu- 
lations on a very dense /c-mesh using the full-potential 
nonorthogonal local-orbital minimum-basis scheme [ p^ . 
Thus we get more detailed information on the Fermi ve- 
locities compared with published data . 

The Fermi surface of MgB2 (see Fig. ^Jconsists of four 
disjoint sheets: two hexagonal heavy hole tubes around 
the T — A line and a light hole and an electron honey- 
combs. In the spirit of a two-band model we divide the 
quasiparticles into subgroups with approximately similar 
material parameters. We choose the heavy-hole tube as 
the first effective band (central blue piece in the middle 
of Fig. 0). Its DOS contribution iVh(O) = iVi(O) is about 
a third of the total DOS: A^(0)/iVi(0) «3 and we adopt 
iV2(0)/A^i(0)=2. The average Fermi velocity over this 
hole tube is Wi?i=4.5x lO'' cm/s while for the remaining 
fast Fermi surface sheets we choose Wi?2 = llxlO^ cm/s. 
In the clean limit the scattering rates ^imp,ij can be set 
to 0. Lets start with weak-coupling scenario. If one 
could ignore the fast quasiparticles entirely, one would 



return to the ISB case, considered above, with A = 0.8 
and ?;f=4.5x10'^. From Eq.l one then gets iJc2(0)=6.5 T 
which is still much smaller than the experimental value. 
Moreover, within a consequent isotropic coupling sce- 
nario the interband coupling cannot be ignored. As a 
result, we arrive at 2T, only (see Fig. 2). For compari- 
son the ISB case for the same A and vp~^.9 x lO^cm/s 
is shown, too. For this moderate anisotropy of vp, the 
two-band effect is practically unimportant. Hence, to 
describe Hc2{T) quantitatively, an enhanced coupling of 
the heavy holes would be needed. Deformation potential 
calculations provide support for such an enhancement for 
the larger heavy- hole tube [pl| . 

It was shown above, that at given Tc = 40 K and el-ph 
interaction with high-frequency phonons, only, the cou- 
pling constant cannot exceed unity and i?c2(0) doesn't 
exceed several Tesla at most. To increase the cou- 
pling constant, we are forced to adopt strong coupling 
to low-frequency mode(s). Hence, we adopt a two-peak 
spectral function containing for instance a low-frequency 
peak at 17.5 meV as well as an effective high-frequency 
peak at 60 meV with comparable spectral weights. We 
adopt Agomev/AiT.smey^O.SS. Then reasonable agree- 
ment with available experimental data for Hc2{T) [p^jl7 
can be achieved with Aii=1.5, A22=0.4 and A2i=0.5 (see 
Fig.2). With the bare band value N(0)=1.7 mJ/mol 
]ll| we arrive at 75=4.1 mJ/mol and ujpi=7.5 eV 
in accord with experimental data |l0|jT^Jl^ . The calcu- 
lated gap values 2Ai=5A:bTc and 2A2—3kBTc above and 
below the BCS value S.bksTc, respectively, can be com- 
pared with the experimental values taken from tunneling 
(SAkBTc, iMkBTa) (l9[|2|l and specific heat {2AkBTc) 
ll and NMR data {2.5kBTc,5.kBTc) 

It is interesting to realize, that switching off the soft 
mode at fixed remaining parameters results in Tc^hard ~ 
7K, only. However, switching off instead the high- 
frequency modes, we would arrive at Tc.soft ~ 20K at 
least. Alternatively, if the spectral functions differ signif- 
icantly for various Fermi surface sheets, the contribution 
of the high-frequency modes to Tc can enhanced up to 
about 20 K. Anyhow, this means, that the "medium" 
high Tc superconductivity of MgB2 within our model 
is not provided by a simple "hydrogen-like" scenario as 
frequently discussed in the literature, but by a more 
complex scenario, where both Tc and i?c2(0) do depend 
very sensitively on the coupling with the soft-mode. A 
somewhat similar interplay of high- and low-frequency 
phonons occurs in the well-known PdH(D) palladium hy- 
drides (deuterides). 

The strong coupling in the heavy-hole band due to in- 
teraction with the soft mode needs some interpretation 
and discussion. A simple down-tuning of the calculated 
i?2g-mode frequency from 58.5 meV [|l2| to 17.5 meV, pre- 
serving thereby the harmonic approximation, would re- 
sult in large zero-point fluctuations amplitudes for which 
there is not enough space in a real crystal. In addition 
we would arrive at an unrealistic super strong coupling 
limit: Ai = 10.4 adopting the estimated Hopfield parame- 



ter 771 = 8.2eV/A^ These difficulties can be resolved 
adopting an anharmonic model for the lattice vibrations. 
Within its extreme anharmonic two-level system approx- 
imation (deep double- or multiple wells) ||2^,Q the cor- 
responding anharmonic electron-boson coupling constant 
reads 
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Ai(r) — —riidQta.nh{hLLji/2kBT)/huJi, 



(4) 



where do ~ 0.1 A is a typical double- well distance. It 
is sizable and at the same time not too large: thus we 
arrive at Ai=2 in accord with An + A12 obtained above 
phenomenologically. Quite interestingly, in this limit the 
boson frequency remains constant whereas the intensity 
of its spectral density increases with decreasing T like 
in the experimental neutron data |^ until saturation at 
low T is achieved. For the adopted value of uJi—17.5 
meV this occurs below 50 K, i.e. only slightly above Tc- 
Alternatively, in principle, a softening of other modes 
than the i?2g mode should be considered, too. Guided 
by the experience obtained from the somewhat related 
transition metal borocar bides (see below), a softening of 
the acoustic modes ||2^ is worthy to be considered. On 
the other hand, despite the soft mode the general features 
of the phonon density of states till 60 meV are reasonably 
well described by the model JTsf . 

There is an interesting similarity between the picture 
we arrived at above for MgB2 and the picture based on 
the present knowledge on transition metal borocarbides 
such as YNi2B2C. The latter exhibit pronounced nesting 
properties of a rather complex Fermi surface (FS). There 
the strongly coupled slow electrons stem from the nested 
parts of the FS. The anisotropy of the FS in borocar- 
bides is somewhat larger than in the present case but the 
anisotropy of the coupling in MgB2 seems to be larger. 
In this context the g-dependent measurement of the soft 
phonon near 17 meV in the superconducting state is of 
great interest. Analogous behavior as for the anomalous 
7 meV mode for YNi2B2C [^sj might be expected. The 
ffc2(0), the positive curvature of Hc2{T) near Tc, and the 
sizeable deviations of the shape of Hc2 {T) at low T from 
the parabolic WHH (Werthamer-Helfand-Hohenberg) - 
shape in both systems can be ascribed to the interband 
coupling with the remaining weakly coupled fast quasi- 
particles. 

In conclusion, the proposed ISB criterion excludes the 
applicability of the frequently used isotropic single band 
model. The shape and the magnitude of the upper crit- 
ical field in MgB2 can be successfully described within 
a mute-band Eliashberg model with two options. Pro- 
vided the modern band structure calculations describe 
correctly the electronic structure near the Fermi surface, 
then there should be a sizeable coupling to some soft 
bosonic modes, not necessarily phonon-like, not predicted 
by harmonic phonon calculations but seen possibly in the 
neutron data. Otherwise the bare quasiparticle Fermi 
velocities should be considerably smaller than the band 
structure predictions. 
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TABLE I. ISB check for selected superconductors. Input 
quantities: ujpi in eV, i/c2(0) in Tesla, cell volume V in A'^, 
Sommerfeld constant 75 in mj/mol K^, Tc in K, coupling 
constant A, and the ISB parameter Q (Eq. (1)). Row 3: weak 
coupling, (Q) —60 meV; row 4: strong coupling (r2)=17.5 
meV 
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"Ref. 26[|27| and references therein. 
''Ref. [ 7|] and references therein. 

"Refs. ]l^ , p^ , p^ , weak coupling scenario, (r2)=60 meV 
''ibid, strong coupling scenario, (r2) = 17 meV 
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FIG. 1. (Color) The Fermi surface of MgB2 in the full 
Brillouin-zone. The P-point is the midpoint of the hexag- 
onal prisma and the A-point(s) are the midpoints of the 
hexagons at the top and the bottom of the prisma, respec- 
tively. Light-hole tube (upper panel); Heavy-hole tube in- 
side the hole honeycomb (middle panel); Electron honeycomb 
(lower panel). The Fermi velocities are measured in units of 
10^ cm/s presented by different colors shown at the bottom 
of the figure. 



□ exp.data Bud'ko et al 




FIG. 2. Experimental data for Hc2{T) of MgB2 com- 
pared with the theoretical curve of the two-band model for 
both anisotropic coupling and Fermi velocities (full line). Pa- 
rameter set see text. For comparison the case of a com- 
plete isotropic single band model (ISB) (dashed) and the 
two-VF-band model with isotropic coupling (dotted line) are 
shown, too. 



